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The collapsar model is the most promising scenario to explain the huge release of energy 
associated with long duration gamma-ray-bursts (GRBs). Within this scenario GRBs are 
believed to be powered by accretion through a rotationally support torus or by fast rotation 
of a compact object. In both cases then, rotation of the progenitor star is one of the key 
properties because it must be high enough for the torus to form, the compact object to 
rotate very fast, or both. Here, we check what rotational properties a progenitor star must 
have in order to sustain torus accretion over relatively long activity periods as observed in 
most GRBs. We show that simple, often cited, estimates of the total mass available for torus 
formation and consequently the duration of a GRB are only upper limits. We revise these 
estimates by taking into account the long term effect that as the compact object accretes 
the minimum specific angular momentum needed for torus formation increases. This in turn 
leads to a smaller fraction of the stellar envelope that can form a torus. We demostrate that 
this effect can lead to a significant, an order of magnidute, reduction of the total energy 
and overall duration of a GRB event. This of course can be mitigated by assuming that the 

' 
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progenitor star rotates faster then we assumed. However, our assumed rotation is already 
high compared to observational and theoretical constraints. We estimate the GRB duration 
times, first by assuming a constant accretion rate, as well as by explicitly calculating the 



$_i 1 free fall time of the gas during the collapse. We discuss the implications of our results. 



1. Introduction 

The collapsar model for a gamma ray burst invokes a presence of an accretion torus 
around a newly born black hole (Woosley 1995; Paczyhski 1998; MacFadyen & Woosley 
1999). The accretion energy is being transferred to the jet that propagates through the 
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collapsar envelope and at some distance from the central engine is responsible for producing 
gamma rays. This type of model is commonly accepted as a mechanism for a long gamma 
ray burst production, because the whole event can last as long as the fallback material from 
the collapsar envelope is available to fuel the accretion disk or torus. 

However, one should bear in mind that the rotating torus may form only when the 
substantial amount of specific angular momentum is carried in the material (see e.g. Lee & 
Ramirez- Ruiz 2006 for a recent study of this problem). This can be parameterized by the so 
called critical specific angular momentum value, which is dependent on the mass of the black 
hole, i.e. Z cr i t = 2R g c, where R g is the gravitational radius. Because the black hole mass is 
not constant during the collapsar evolution, but increases as the envelope material accretes 
onto it, the critical angular momentum will change in time. Consequently, the amount of 
the rotating material, which was initially available for the torus formation, may become 
insufficient at a later stage of the collapsar evolution. 

Moreover, the spin of the black hole will be changed by accretion. Whether the black 
hole can achieve a spin parameter close to the maximum one, depends on the properties of 
the accreted mass. While a large spin (a ~ 0.9) is thought to be a necessary condition for the 
jet launching (Blandford & Znajek 1977), it may happen that not enough specific angular 
momentum is being transferred to the black hole as its mass increases. 

Another challenge for the collapsar model is due to the effects of stellar wind, which 
removes from the Wolf-Rayet stars a large fraction of angular momentum (Langer 1998). 
However, the winds of massive stars are relatively weaker in the low metallicity environment 
(Abbott et al. 1982), and possibly the GRB progenitor stars can rotate faster than an 
average W-R star (Vink 2007). 

Here we address the question of whether the collapsing star envelope contains enough 
specific angular momentum in order to support the formation of the torus. Furthermore, it 
will be interesting to study the problem of spin up and spin down of the newly born black hole, 
and we shall consider this in a follow up paper. These two are the key properties needed 
to launch the GRB jet for an extended period of time. Because the angular momentum 
distribution in the Wolf-Rayet stars is unknown, we may pose this question also in a different 
way: we want to check how much angular momentum has to be initially present in the stellar 
envelope in order to allow for a long GRB production. 

In Section [2], we describe the model of the initail conditions and evolution of the collaps- 
ing star, adopting various prescriptions for the angular momentum distribution and different 
scenarios of the accretion process. In Section [3l we present the results for the mass accreted 
onto the black hole in total and through the torus. We also estimate the resulting GRB du- 
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rations, first in case of a constant m and then by explicitly calculating the free fall velocity of 
the gas accreting in the torus. Finally in Section HI we discuss the resulting duration time of 
GRBs, as a function of the distribution of the specific angular momentum in the progenitor 
star. 



2. Model 

In the initial conditions, we use the spherically symmetric model of the 25 M & pre- 
supernova (Woosley & Weaver 1995). The same model was used by Proga et al. (2003) 
in their MHD simulation of the collapsar. Figure [T] shows the density profile and the mass 
enclosed inside a given radius. The Figure also shows the free fall timescale onto the enclosed 
mass, corresponding to the radius. 

The angular momentum within a star or rotating torus may depend on radius (see e.g. 
Woosley 1995, Jaroszyhski 1996, Daigne & Mochkovitch 1997 for various prescriptions). Here 
we parameterize this distribution to be either a function of the polar angle 9 (simplest case; 
models A and B), or a function of both radius r and 9 (models C and D). 

First, we assume the specific angular momentum to depend only on the polar angle: 

4pec = kf(9). (1) 

We constitute two different functions: 

f(9) = 1- |cos0| (model A) (2) 

f(0) = sin 2 9 (model B) (3) 
The rotation velocity is therefore given by: 

<V = —J{0) (4) 
r sin 9 

The normalization of this dependence is defined with respect to the critical specific angular 
momentum for the seed black hole: 

l = x/ crit (M° H ) = x x 3.54 x io 16 M^®l cm 2 s -i ( 5 ) 

where R g = 2GMg H /c 2 is the Schwarzschild radius (non-rotating black hole). 

Second, we assume that the specific angular momentum will depend on the polar angle, 
as well as on the radius in the envelope, as: 

4 pec = lo9{r)f(9), (6) 
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Fig. 1. — The density and mass profiles in the pre-supernova model. The data are taken 
from Woosley & Weaver (1995), model No. S251S7B@14233. The x-axis on the upper panel 
shows the free fall timescale corresponding to the radius shown on the lower panel's x-axis. 
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We adopt the following functions: 

4 P ec=^crit( ) 2 sin 2 6» (model C) (7) 

'"core 

^spec = x \/SGM~r sin 2 6 (model D) (8) 

The above model C corresponds to the rotation with a constant angular velocity fl, while the 
model D corresponds to a constant ratio between the centrifugal and gravitational forces. 
Note that the strong increase of / spec with radius will lead to a very fast rotation at large 
radii. Therefore, a cut off may be required at some maximum value, / max (see below). 

The normalization of all the models is chosen such that the specific angular momentum 
is always equal to the critical value at 9 = 90°, and at r = r core if the model depends on 
radius. In Section [3], we present the results of our calculations considering a range of initial 
values of x. 

Initially, the mass of the black hole is given by the mass of the iron core of the star: 

/"^corc 

M° BH = M corc = 4tt / pr 2 dr. (9) 
J o 

For a given x, a certain fraction mass of the collapsar envelope, Mf, carries a specific angular 
momentum smaller than critical Z° rit = Z C rit(^BH") : 

M° = 2tt / Pl r 2 sm6d6dr (10) 

J f core " 

where p\ = p(r, 9) | M o , is the density in the envelope where the specific angular momentum is 
smaller than critical. Here, the radius r max is the size of the star. Correspondingly, by M 2 ° we 
denote the fraction of the envelope mass that carries the specific angular momentum larger or 
equal to the critical, with p 2 = p(r, 0)\ t>l o . , and the total envelope mass is M e ° nv = M° + M 2 ° . 
Only the mass M° can form the torus around the black hole of the mass Mg H . 

The above relations set up the initial conditions for the torus formation in the collapsar, 
and l crit is defined by the mass of the iron core, M core . However, as the collapse proceeds, 
the mass of the black hole will increase and the critical specific angular momentum will be a 
function of the increasing mass: / C rit(^fBH)- The main point of this work is to compute the 
mass of the progenitor with / > / crit , taking into account this effect. Below, we redefine pi 
and P2, so that the Z cr it — A^br" relation is taken into account. 

To compute the mass of the envelope part that has high enough / to form a torus around 
a given BH, and to estimate the time duration of the GRB powered by accretion, we need 
to know the mass of this black hole. A current Mbh depends on the mass of the mass of 
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the seed black hole and the accretion scenario. We approximate this scenario in a following 
way. We assume that accretion is nonhomologous and BH grows by accreting mass Am k , 
which is a function of the mass of a shell between the radius ?\ and r\ + A?\ (see e.g. Lin 
& Pringle 1990). Formally, we perform the calculations of Mrh and Am k iteratively: 



Here p depends on the accretion scenario (see below) and contains the information of the 
specific angular momentum distribution. The above two equations define an iterative pro- 
cedure due to the nonlinear dependence of p on M B h- We start from the radius r = r core , 
i.e. that of an iron core. 

We distinguish here three possible accretion scenarios: 

(a) the accretion onto black hole proceeds at the same rate both from the torus and from 
the gas close to the poles, with I < Z k rit , i.e. p = p (and does not depend on 9); 

(b) the envelope material with I < Z k rit falls on the black hole first. Thus, until the polar 
funnel is evacuated completely, only this gas contributes to the black hole mass, i.e. p = p\. 
After that, the material with I > Z k rit accretes, and p = p 2 ; 

(c) the accretion proceeds only through the torus, and only this material contributes to the 
black hole growth i.e. p = p 2 . In this case the rest of the envelope material is kept aside 
until the torus is accreted. 

The densities p\ and p 2 , defined above, depend on / k rit = I cr i t (Mg H ). The above accretion 
scenarios are illustrated in the Figure [2j The panel (a) shows the scenario of a uniform 
accretion, in which the whole envelope material falls into black hole, regardless of its specific 
angular momentum. The red color marks the material with I < Z k rit . The blue colors mark 
the material with / > Z k rit , and when the black hole is small, this condition is satisfied for 
a larger range of 9 (dark blue). When the black hole increases, the material with I > Z k rit 
occupies narrower 9 range (light blue). The panel (b) shows the scenario with two steps: 
first the material with I < Z k rit accretes onto the black hole, increasing its mass; after this 
material is exhausted, the material with / > Z k rit starts accreting. Because the black hole 
mass has already increased, material with large / is concentrated very close to the equator. 
The panel (c) shows the scenario in which only the material with / > Z k rit accretes. 

In scenario a the mass accretion rate does not depend on the specific angular momentum. 
This is a very special and simple scenario. In reality, the accretion rate can depend on the 
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where the increment of mass of the black hole is : 




(12) 
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specific angular momentum. For example, if an accreting torus produces a very powerful 
outflow, the weakly rotating polar material could be expelled and never reach the black 
hole (scenario c). This would also be a special and extreme situation. It is more likely 
that both the polar and disk material accrete but at different rates. However, it is unclear 
what these rates are and detailed simulations of MHD flows show that the rates can depend 
on time. For example, there are periods of time when the polar material accretes faster 
than the disk material and vice versa (e.g. Proga & Begelman 2003). To bracket this 
more realistic situation, we consider here another extreme and rather artificial scenario b, 
which corresponds to an somewhat 'reversed' scenario c. In this scenario, initially the torus 
accretion rate is zero and accretion is dominated by the polar material. Only after the polar 
mateial is exhausted, the torus accretion starts. We note that although this scenario is quite 
extreme, it may be relevant if jets in GRBs must be very clean and light because in this 
scenario jets will be moving in the 'empty' polar funnels. 

Due to the increasing mass of the black hole, the critical angular momentum also in- 
creases, and as a result less and less material can satisfy the condition for the torus formation 
(/ > /c r i t ). We stop the calculations, when there is no material with I > l^ Tit , i.e. able to form 
the torus: 

M k 27r f max fJ p 2 r 2 sin OdOdr 

Wk = J^= Jr ^ J r y 2 = o. (is) 

Alternatively, the iterations may be stopped earlier, for example if we impose a physical limit 
based on the free fall timescale or the accretion rate, to be adequate to power the prompt 
GRB phase. 

The duration of the GRB could be estimated as the ratio between the mass accreted 
through the torus, and the accretion rate fa: 

h 

"-max 

M^ir = J2 M 2 ( 14 ) 

k=l 
A/rtoms 

At GRB = (15) 
m 



where the number k max is defined by the Equation 

Note that we assume here the GRB prompt emission is equal to the duration of the torus 
replenishment. In principle, rh may depend on time. Here we take two approaches. First, 
for simplicity, we assume a constant accretion rate of a moderate value (m = 0.01 — 1.0 M Q 
s^ 1 , see e.g. Popham, Woosley & Fryer 1999; Janiuk et al. 2004). Second, in more detailed 
calculations we determine the instantaneous accretion rate during the iterations, determined 
by the free fall velocity of gas in the torus. 
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3. Results 

3.1. Models with the specific angular momentum dependent only on 9 

The Figure [3] shows the initial fraction of the envelope mass which contains large enough 
angular momentum to form the rotating torus, Wq = M^/M^ (see Eq. ITBl for models A 
and B. For instance, for the adopted function f(9) given by Eq. [2] (model A), and for the 
initial angular momentum normalization of x = 1.15, we obtain wq = 0.13. This means that 
only 13% of the total mass of the envelope will initially be able to form the torus, while the 
remaining 87% of the envelope mass will fall radially towards black hole. On the other hand, 
for x > 5, more than 75% of the envelope mass will be rotating fast enough to contribute to 
the torus formation. The model B gives systematically larger values of Wq and for x = 1.15, 
we have w = 0.36, while for x > 5 we have more than 85% of the envelope mass able to 
form the torus. 

As we show below, these are only the upper limits for the mass that could be accreted 
through the torus, and drive the GRB duration. These values will be much smaller, when 
we calculate the collapsar evolution with increasing Z^ rit instead of Z° rit . 

The Figure 0] shows /„ it , i.e. Z crrt as a function of the current radius r^, which is the 
inner radius of the collapsing envelope in a subsequent step k. The figure thus shows how the 
critical specific angular momentum changes with time during the collapse, for an exemplary 
value of x = 7. 

The / cr it rises with time, as the black hole accretes mass from the envelope, and cor- 
responds to a changing black hole mass. The most steep rise is for the uniform accretion 
scenario a, and in this case by definition the result does not depend on the adopted distri- 
bution function for the specific angular momentum, f(9). Therefore both curves marked by 
a solid line overlap. Also, in scenario a the plotted curves do not depend on x, as well as 
neither on the slope nor on the location of the curve. The latter influences only the max- 
imum of this curve, as for larger x we have more material available to form the torus. In 
particular, for x = 7, the two overlapping curves shown in the figure end at r ~ 10 13 cm. 

For the scenario c, i.e. accretion of gas with I > l^ rit , / crit rises less steeply with r k than 
in scenario a, because now less material is contributing to the black hole mass. In this case 
f(9) affects the results, and model A gives systematically smaller values of Z cr it than model 
B. 

For the scenario b, the result is very sensitive to x, and we can have either one or two 
phases of accretion: only the polar inflow or first the polar inflow and then torus accretion. 
The value of x = 7 was chosen, because in model A still no torus is able to form, and we 
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have only phase 1, while in model B this value of x is already large enough and the phase 2 
occurs. 

For phase 1 in scenario b (marked by the thinner lines in the figure), i.e. the material 
with I < l^ rit is accreting, the dependence on f{9) is the following: model A adds more mass 
to the black hole and therefore it leads to the larger values of Z cr it than model B. For phase 
2 of scenario b (present only for model B and marked by the thick line in the figure), the 
evolution starts from the last l^ rit achieved in the end of phase 1. Then l^ rit increases, and 
ultimately reaches the final solution of models Be and A b, because this l cr \ t corresponds 
to the black hole mass that has increased maximally: either only through a torus, or first 
through the through the polar funnels and then through the torus accretion. 

All the curves in Figure H] exhibit a characteristic evolution of their slopes, tracing the 
density distribution in the progenitor star (see the top panel of Figured]). First, the fast rise 
is due to accretion of the most dense inner shells of the stellar envelope. Then, the slope 
flattens, as the density in the envelope decreases and the mass does not grow very fast. In 
the end, the slope of Z„ it = / CT it( r k) = Z cr i t (M B H) rises again, due to larger volume of the 
shells, but this rise is depending on the adopted scenario. In scenario c the sequence is the 
following: increase of Z^ rit — > more accretion — > larger increase of Z^ rit less accretion. In 
the phase 1 of scenario b, such an equilibrium is not established, because the accretion onto 
black hole proceeds through the polar funnels, i.e. using the material with I < Z^ rit , so we 
have: increase of Z^ rit — > more accretion — > further increase of Z~ rit , and for large radii the 
slope of the curves shown in Fig. H] in this scenario is much steeper. The phase 2 of scenario 
b results in the changes of the slope of Z^ rit similar to other scenarios, a and b (provided that 
the phase 2. occurs). 

In Figure we show the total mass accreted onto a black hole, and in Figure we 
show the mass accreted through the torus, both as a function of x. Figure can be also 
regarded as showing the estimated duration of a GRB, if the accretion rate is constant. The 
results are again presented for the 3 scenarios of accretion: a, b and c (marked by solid, 
short-dashed and long dashed lines, respectively), as well as for the two prescriptions for the 
function f(6), models A and B (marked by squares and circles, respectively). The upper 
panels in Figs. and [H] show the results obtained is case of the maximum limit for the free 
fall time set yo t^ ax = 1000 s. 

The values of the total accreted mass are the largest for the scenario of the uniform 
accretion, a. Depending on x, the black hole mass is growing until there is no further 
material with I > l^ lit . 

In scenario b, initially we add to the black hole mass (thus increasing Zjf rit ) , only the 
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material with I < £jf rit . For small values of x, the total accreted mass is the same as for 
scenario a, because the process of accretion lasts in both cases until M\ = 0, i.e. the 
envelope contains no further material with I > l^ Tit . However, for large x (> 7 in model A 
and > 5 in model B), after the black hole has swallowed the whole funnel with I < l^ Tit , 
there is still some material with large specific angular momentum and phase 2 occurs. The 
material accretes now through the torus, but only as long as it has / > l^ Tit . Therefore, for 
large x, the total mass which the black hole gains in scenario b is less than in the scenario a. 

Scenario c assumes, that the BH accretes only the material with I > l^ Tit . Now, the total 
accreted mass can be either the same (for small x) or smaller than in scenario a. 

In Figure M, we show the accreted mass which had I > Z„ it . This represents the accretion 
through the torus, and may be regarded as a direct measure of the GRB duration. Scenario 
a results in a linear scaling of M*°™ s with x: 

M^ s = ax + (3 (16) 

and with a linear fit we obtained a ~ 0.83, (3 ~ —1.41 in model A, and a ~ 1.12, (3 ~ —1.55 
in model B. 

Scenario b predicts that the torus accretion is possible only for large x, while for small 
x torus will not form, as discussed above. The scenario c, by definition, predicts the torus 
accretion for any value of x > 1.0. Therefore, the amount of material accreted with / > l^ Tit 
is in this scenario larger than in scenario a, because the black hole mass grows more slowly. 
Both scenarios a and b result, for large i, in a nonlinear dependence of M*°™ s on x. 

To sum up, in the models A and B, i.e. if the specific angular momentum in the 
progenitor star depends only on the polar angle, the total mass of material capable of forming 
the torus can only be a small fraction of the envelope mass. Depending on the accretion 
scenario, it is at most ~ 3.5M , i.e. 15% of the envelope mass, for l = 3/° rit = 10 17 cm 2 s _1 , 
and between ~ 7 and ~ 15M , i.e. 30%-65% of the envelope mass, for / = 10/° rit = 3.3 x 10 17 

2 —1 

cm s . 

Note that we could proceed with larger / , but we stopped our calculations at x = 10, 
because larger x would already imply very fast rotation at the equator. In the present section 
we did not assume any maximum limit on the specific angular momentum; this will be taken 
into account in the next section. Howewer, we considered the effects of the maximum limit 
for the free fall timescale. As shown in the Figures, the limit of £™ ax = 1000 s plays an 
important role when x > 5, and in all the models and scenarios the dependence of the 
accreted mass on x is significantly weaker than for the case of no limiting tg. For large x, in 
scenario a the total mass accreted onto BH is constant with x and equal to the fraction of 
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the envelope mass enclosed within the radius r w 1.58 x 10 11 (see Fig. [I]). The mass accreted 
via torus is smaller than that, and for x — 10 it reaches about 6 solar masses. 

3.2. Models with the specific angular momentum dependent on r and 9 

Now, we investigate how the total accreted mass and in consequence the duration of the 
GRB will be affected if the angular velocity Q in the collapsing star is constant or given by 
a fixed ratio of the centrifugal to gravitational force (Equations [7] and [H respectively). 

In these two models, C and D, the specific angular momentum is a strong function 
of radius. Therefore, if we do not impose any limit on the specific angular momentum, 
Iraax, the outer shells of the star will always have a substantial amount of specific angular 
momentum, larger than any current critical value. Consequently, the GRB will continue 
until the last shell is accreted. However, this would imply a very fast rotation of the star in 
its outer layers. This would be inconsistent with the progenitor models for GRBs, so in a 
more realistic version of our modeling we impose a maximum value of the specific angular 
momentum. This assumption will lead to a shorter GRB event, as it will end when the 
increasing black hole mass will imply i^St > ^max- 

First, we investigate how the black hole mass, Mg H , and critical specific angular mo- 
mentum, /k rit , depend on the accretion scenario. For the scenario a (i.e. uniform accretion), 
the black hole is fed by the whole envelope regardless of the local specific angular momentum 
value. The result is the same as in the cases explored in Section 13.11 the total accreted mass 
does not depend neither on the distribution function (C or D) nor on the normalization (x). 
In the Figure we show how the critical specific angular momentum increases when the 
subsequent envelope shells are accreting (i.e. as a function of radius). The solid lines in this 
figure overlap, and basically the curve is the same as in the figure H] (for models A and B), 
the only difference being that the maximum value reached in models C and D can be larger 
(specifically, for log/ max > 17.3). 

The amount of the total accreted mass in this case is constant (Figure [8]) and the value 
depends only on Z max : 

M accr = ^ - M core . (17) 

For our model of the star and Z max = 10 17 cm 2 s _1 this gives M accr = 3.99M Q . If there is no 
cutoff of Z max then simply the total envelope mass is accreted, 23.9 M Q (cf. the bottom panel 
of Fig. ED- 

The situation becomes more complicated when we adopt the scenario c (accretion of 
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material with I > i^ t ). The accreted mass in this scenario depends both on the distribution 
and on the normalization of the specific angular momentum in the pre collapse star. 

For small x, the accreted mass is small, because the process ends when I < l^ rit every- 
where in the envelope and M\ = 0. The total mass accreted on black hole is sensitive to 
the model distribution function. In particular, the fact that this function strongly depends 
on the radius means that the inner shells contain mostly material with I <C ^ rit . Thus, only 
the more distant shells contribute to the mass of the black hole (see Fig. [7J dashed lines), 
and the particular value of for which the black hole mass and /^ rit start rising depends on 
x. Note that Fig. [TJ shows only the results for x = 0.05 (arbitrarily chosen value). 

For large x, the accreted mass will asymptotically reach the result for the scenario a 
(see Figure El bottom panel), because if x — 1 then the whole envelope material satisfies 
I > ^crit- Therefore the Z C rit(^BH") functions seen in the Figure [7] will eventually overlap if x 
is close to 1.0. Therefore in the Figure [8] we show only the results for x < 1.0 because these 
are the most interesting: for larger x the mass accreted both through the torus and in total, 
approaches a constant value. 

The smallest amount of the total accreted mass is obtained when we impose a cutoff 
limit on the specific angular momentum, Z max . This is shown in Fig. [8] (middle panel) for 
the value of Z max = 10 17 cm 2 s -1 . The total accreted mass is in this model M*°* r <C M env , and 
very weakly depends on x. 

The value of the / max cut off has to be chosen carefully, because if / max > ^-M env = 
4.23 x 10 17 cm 2 s~ 1 , then the accreted mass would be equal to the envelope mass (and equal 
to that obtained with the uniform accretion scenario), for x = 1. Any value of Z max larger 
than the above value will not affect the results. 

The chosen form of the specific angular momentum distribution (models C or D) only 
slightly affect the results. For x ~ 0.01, the model C gives larger value of accreted mass, 
while for x > 0.1, the model D leads to somewhat larger M*°*. However, the results in this 
case are also affected by the numerical issues because of the very small number of steps after 
which the calculation is finished. 

The accreted mass will be zero, and the burst will not occur, if the normalization x is 
very small. The minimum value can be estimated for C and D models separately, if we take 
the specific angular momentum to be everywhere smaller than critical: 



,c 



r, 



core \ 2 



) 2 = 1.22 x 10 



-li 



(18) 



x. 



mm 



r, 



ma 



and 




(19) 



13 



Now, we can estimate the duration of GRB by means of the mass accreted onto the black 
hole with I < l^ rit , i.e. via the torus. In scenario c, this is, by definition, the same as the 
total mass accreted. As can be seen from the bottom panel in the Figure EE M*°™ s > 2OM , 
and GRB duration on the order of ~ 40 seconds, is possible only with the model with no cut 
off on / max and x > 0.1. For both models C and D the uniform accretion scenario a gives 
slightly less amount of mass accreted through the torus, but the difference is visible only for 
x < 0.1. 

The more physical models with the angular cut off to Z max = 10 17 cm 2 s _1 always give 
less than ~ 4M Q of mass accreted via torus, which corresponds to the GRB duration of only 
about 8 seconds. In scenario a, the mass accreted with I > Z^ rit is even smaller than that, 
especially for small x (< 0.5). No mass will be accreted through the torus if x < 0.05 (model 
C) or x < 0.1 (model D). 

In scenario c, the mass accreted with I > l^ Tit is independent on x, if there is a cut off 
on Z max (cf. Eq. [IT] and note that this mass is equivalent to the total mass accreted). 

The accretion scenario b, i.e. that of the accretion composed from two steps: polar 
funnel and than torus, is not discussed for models C and D. It is because only a very small 
fraction of the envelope, and for very small x, has I < l CT i t , so basically the results will not 
differ much from the scenario c. 

Finally, we tested the models C and D with an upper limit imposed on the free fall 
timescale, tg ax = 1000 s. To comaper the two effects, in these tests we did not assume any 
limit for the specific angular momentum / max . As shown in the upper panels in the Figures 
[8] and [9j the mass accreted on the black hole (both in total and via torus) is now 3 times 
smaller than without the t^ ax limit. Nevertheless, this effects is not as strong as the limit 
for the progenitor rotation law imposed by / max . 

To sum up, in the models C and D, i.e. if the specific angular momentum in the 
pre-collapse star depends on 9 and r in such a way, that either the angular velocity Q is 
constant, or constant is the ratio between the gravitational and centrifugal forces, a fraction 
of envelope material able to form a torus can be much smaller, or much larger than in models 
A and B. The fraction of 100% is possible if there is no limiting value on the specific angular 
momentum (or, more specifically, the limiting value exceeds 4.23 x 10 17 cm 2 s _1 in our model), 
and no limit for a free fall timescale. However, in more physical modeling which accounts 
for such limits, this fraction becomes very small: for tg ax = 1000 s we obtain about 30% of 
the envelope accreted via torus, and for Z max = 10 17 cm 2 s _1 we obtain at most 15%. 
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3.3. Duration of a GRB 

In Figure [10] we show the instantaneous accretion rate during the iterations, i.e. as 
a function of the current radius. As the Figure shows, the accretion rate is the largest at 
the beginning of the collapse, and equal about 0.1 M & s _1 . In model C, for x = 0.05 the 
condition for torus formation is not satisfied initially (cf. Fig. [7]), so the accretion rate 
through the torus is zero. For the same x, in model D the torus is already formed near the 
equatorial plane, and the accretion rate is about 0.03 M s _1 . 

The duration of a gamma ray burst depends on the accretion rate and the total mass 
accreted via torus onto the black hole. If the torus is a poer source for a GRB, the accretion 
rate, albeit must not be constant, cannot drop to a very small value either. Below, say 0.01 
solar masses per second the neutrino cooling in the accretiong torus may become inefficient. 
Table [1] shows the results for the computations, in which we have limited the iterations to 
the minimum accretion rate of rh min = 0.01M Q s _1 . The table summarizes the mass acreted 
via torus for all of our models for the progeniotor rotation (i.e. A, B, C and D) as well as 
the two accretion scenarios (a and c). The mass accreted through the torus never exceeds 4.5 
M & s _1 , thus implying that the limiting accretion rate value influences the results stronger 
than the limits for the free-fall time, and comparably to the limit for progenitor rotation. 

The Table [2] summarizes the durations of GRB prompt phase, also for the four models 
and two accretion scenarios. The results correspond to the total masses accreted via torus 
that are given in Table [TJ The duration time was calculated as tcRB = M^°l us / <fh>, where 
< rh > is the mean accretion rate during an iteration. Note, that because the minimum 
accretion rate was fixed at 0.01 solar masses per second, the average value is equal to 0.5rh max . 

The calculated duration times of GRBs are the largest for models C, because the average 
accretion rate in these models is the smallest. Taking into account only the free fall timescale 
and under the adopted assumptions for a minimum rh, these models give at most ~ 145 s of 
the GRB prompt phase. All the other models result in the GRB duration below 50 seconds. 

4. Discussion and conclusions 

The durations of GRBs range from less than 0.01 to a few hundred seconds (for a review 
see e.g. Piran 2005), and the long duration bursts, T 90 > 2 s, are supposed to originate from 
the collapse of the massive rotating star. The collapsar model assumes that the presence of a 
rotating torus around a newly born black hole is a crucial element of the GRB central engine 
for the whole time of its duration. In this work, we found that some specific properties of 
the progenitor star are important in order to support the existence of a torus, which consists 
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Table 1: The mass accreted through the torus (in solar masses) for various models and accre- 
tion scenarios, under the assumption that the minimum accretion rate is rh m [ n = O.OIMqS" 1 



X 




Mj£ us 






2.0 


0.35 


0.78 


0.78 


1.37 


3.0 


1.01 


1.58 


1.76 


2.37 


4.0 


1.62 


2.12 


2.61 


2.94 


5.0 


2.14 


2.54 


3.09 


3.27 


6.0 


2.53 


2.78 


3.38 


3.48 


7.0 


2.82 


3.02 


3.57 


OCA 

3.64 


8.0 


3.02 


3.19 


3.73 


3.78 


9.0 


3.19 


3.30 


3.81 


3.85 


10.0 


3.33 


3.43 


1 88 
o.oo 


3.91 


X 


m Ca 


71 /[torus 
m Cc 


7\ /[torus 
m Da 


7i /[torus 
Ivl Dc 


0.05 


1.96 


2.55 


1.64 


2.50 


0.1 


2.93 


3.18 


3.41 


3.53 


0.2 


3.55 


3.66 


4.05 


4.07 


0.3 


3.80 


3.85 


4.21 


4.23 


0.4 


3.95 


3.97 


4.35 


4.31 


0.5 


4.04 


4.09 


4.39 


4.35 


0.6 


4.11 


4.15 


4.43 


4.43 


0.7 


4.17 


4.21 


4.45 


4.45 


0.8 


4.22 


4.26 


4.47 


4.47 


0.9 


4.27 


4.30 


4.48 


4.48 


1.0 


4.32 


4.33 


4.49 


4.49 
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Table 2: Duration of GRB (in seconds) for various models and accretion scenarios, under the 
assumption that the minimum accretion rate is m m i n = O.OlM s _1 



X 


t Aa 


tAc 


tBa 


t Bc 


2.0 


7.00 


15.23 


11.17 


19.75 


3.0 


15.15 


23.77 


21.96 


29.60 


4.0 


21.96 


28.86 


30.72 


34.68 


5.0 


27.32 


32.43 


35.25 


37.33 


6.0 


31.01 


34.08 


37.80 


38.93 


7.0 


33.50 


35.99 


39.46 


40.24 


8.0 


35.28 


37.34 


40.91 


41.49 


9.0 


36.66 


37.96 


41.39 


41.85 


10.0 


37.70 


38.85 


41.90 


42.27 


X 


t C a 


t C c 


tDa 


tDc 


0.05 


106.25 


139.75 


87.84 


99.89 


0.1 


126.99 


147.09 


48.21 


49.95 


0.2 


136.21 


150.19 


47.62 


47.70 


0.3 


139.29 


150.11 


47.08 


47.29 


0.4 


141.32 


150.17 


47.41 


46.99 


0.5 


142.43 


151.77 


47.44 


46.98 


0.6 


143.26 


151.66 


47.32 


47.37 


0.7 


144.21 


142.71 


47.38 


47.38 


0.8 


131.66 


126.07 


47.22 


47.20 


0.9 


117.38 


114.70 


47.30 


47.31 


1.0 


108.81 


106.52 


47.26 


47.28 
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of the material with specific angular momentum larger than a critical value I > l CI it. We 
studied how the initial distribution of specific angular momentum inside the stellar envelope 
affects the burst duration, taking into account the increase of the black hole mass during the 
collapse process. 

Following Woosley & Weaver (1995), we considered the model of pre-supernova star that 
predicts the existence of the ~ 1.7M iron core, which forms a black hole, and the ~ 24M Q 
envelope. Therefore in the simplest approach, when the mass available for accretion is the 
total envelope mass, and the accretion rate is a constant value on the order of 0.1 — O.5M s _1 , 
the central engine is able to operate for a time required to power a long GRB (i.e. several 
tens to a couple of hundreds of seconds) . 

However, McFadyen & Woosley (1999) in their collapsar model show that most of the 
initial accretion goes through the rotating torus rather than from the polar funnel. Torus 
formation is possible if material with substantial specific angular momentum is present in 
the envelope, initially and throughout the event as well. In this sense the GRB duration 
estimated in the uniform accretion scenario is only an upper limit. 

In our calculations, this upper limit is achieved only if the specific angular momentum 
distribution in the pre-supernova star is a strong function of radius (i.e. g(r) ~ r 2 or 
g(r) ~ and the inner parts have x = /// C rit( r m) ~ 1.0 (for the initial black hole mass 
we have Z° rit ~ 3 x 10 16 cm 2 s _1 in our model), while the outer parts of the star may have 
an specific angular momentum as large as Z max > 4.23 x 10 17 cm 2 s _1 . Both these conditions 
challenge the progenitor star models: they require either a rigid rotation of the star, or a 
huge ratio of centrifugal to gravitational force. The latter, if we want to keep the value of 
^centr/^graw ~ 0.02 (as taken by Proga et al. 2003), would lead to the mass accreted through 
the torus equal to be a fraction of the envelope mass, and a correspondingly shorter GRB 
duration (i.e. one should take x = 0.05, which implies the the GRB duration in our 'mass 
units' to be 20-21 M , cf. Fig. E). 

Furthermore, the progenitor star models used in MacFadyen & Woosley (1999; see also 
Proga et al. 2003), would rather assume a limiting value for the specific angular momentum. 
In our modeling we followed that work and calculated an exemplary sequence of models with 
l max = 10 17 cm 2 s _1 . The results in this case are not promising for a long GRB production: 
the GRB duration in the accreted mass units does not exceed 4 M , which would give an 
event not longer than a hundred seconds and only if the accretion rate is very small. 

The models with the specific angular momentum distribution depending only on the 
polar angle, 8, also yield not very long GRB durations. In this case, the mass accreted with 
I > Z cr it is of the order of 15 M & if the specific angular momentum in the progenitor star is 
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about ten times the critical one (i.e. x ~ 10), and the accretion proceeds only through the 
torus (the latter finds support in MHD simulations such those performed by Proga et al. 
2003). If the accretion proceeds uniformly through the envelope, the GRB duration drops to 
about 10 M Q for the same value of x. Finally, in the scenario when accretion proceeds first 
through the poles and then through the torus (i.e. scenario b as indicated by HD simulations 
performed by MacFadyen & Woosley 1999), there is no GRB for x < 5 — 7 (depending on the 
shape of /(#)), because all the material is swollen by the black hole during the first stage. 
For large x, the resulting GRB duration is in between of the scenarios a and c. We plan 
to consider other progenitor star models such as those computed by Heger et al. (2005) to 
check how our conclusions depend on specific radial and rotational profiles. 

We also investigated the models in which the mass accreted onto BH was limited by the 
free fall timescale or the minimum accretion rate. In case of of the free fall time limited to 
1000 seconds, the mass accreted onto the black hole is much smaller than the total envelope 
mass, and reached up to 8 M Q but for a very fast rotation of the progenitor star. Finally, 
the explicitly calculated duration times of GRB, obtained due to the released assumption 
of a constant accretion rate, were at most 30-150 seconds, depending on the model of the 
specific angular momentum distribution and accretion scenario. 

The effects of accreting non-rotating or slowly rotating matter on a black hole can 
reduce also capability of powering GRBs through the Blandford-Znajek mechanism. In the 
estimations done by McFadyen & Woosley (1999) a-posteriori, i.e. using the analytical 
models of the accretion disk to extend the calculations down to the event horizon, the 
authors calculated the evolution of the black hole mass and angular momentum. The initial 
dimensionless angular momentum parameter of the iron core is taken in their work to be 
either a ini t = or a init = 0.5. However, the black hole changes its total mass and angular 
momentum as it accretes matter (see Bardeen 1977 for specific formulae). In this way, if the 
specific angular momentum supply is substantial, even starting from a = 0, a Schwarzschild 
black hole, the finite amount of accreted mass makes it possible to obtain a — 1. On the 
other hand, the material with very small specific angular momentum, which is present in a 
collapsing star, will spin down the black hole. 

The effect of the evolution of the black hole spin due to the accretion (spin up) and the 
Blandford-Znajek mechanism (spin down) has been studied in Moderski & Sikora (1996). 
Lee, Brown & Wijers (2002) studied the case of GRBs produced after the progenitor star has 
been spun up in a close binary system due to spiral-in and tidal coupling. Recently, Volonteri 
et al. (2005) and King & Pringle (2006) calculated the spin evolution of supermassive black 
holes due to the random accretion episodes. In our model, the black hole spin evolution 
is not episodic, but a continuous process. The calculations of the BH angular momentum 



-19- 



evolution are left for the future work. Such calculations may possibly show that obtaining 
a large BH spin parameter, a ~ 0.9, is rather difficult when a large fraction of the envelope 
material has / <C Z crit . 
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Fig. 2. — The scheme of accretion scenarios. The red color indicates the material with 
I < lent- The blue colors indicate the material with I > l CI i t : darker for smaller black hole 
mass, and lighter for larger black hole mass. Arrows indicate, which material is accreting 
and contributes to the black hole growth. 
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Fig. 3. — The initial mass fraction of material with the angular momentum I > Z crit , as 
a function of the initial normalization of the specific angular momentum distribution, for 
model A (solid squares) and model B (open circles) of the distribution function f(9). 



-23- 




Fig. 4. — The critical specific angular momentum as a function of radius within which the 
envelope mass collapsed onto BH. The results are shown for one exemplary value of the 
initial normalization parameter, x = Iq/Ictu — 7. The models of the distribution function 
f(9) are: A (solid squares) and B (open circles), and accretion scenarios are: a (solid line), 
b (short dashed line) and c (long dashed line). The thin line for model B6 represents the 
results from the phase 1 of accretion (through the poles), while the thicker line B6 represents 
the results from the phase 2 (torus accretion). Note, that for the model A b there is only a 
thin line, and no phase 2, because for x = 7 the torus does not form. Note also, that the 
solid lines (i.e. models Aa and Ba overlap. 
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Fig. 5. — The total mass accreted onto a black hole, as a function of the initial normalization 
of the specific angular momentum distribution. The models of the distribution function f{9) 
are: A (solid squares) and B (open circles), and accretion scenarios are: a (solid line), b 
(short dashed line) and c (long dashed line). The upper panel shows the case when the mass 
accretion is limited by a maximum free fall time of 1000 seconds, while the lower panel shows 
the results for no limiting tff (cf. Fig. CD). 
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Fig. 6. — The accreted mass with I > Z cr i t , as a function of the initial normalization of 
the specific angular momentum distribution. The models of the distribution function f{9) 
are: A (solid squares) and B (open circles), and accretion scenarios are: a (solid line), b 
(short dashed line) and c (long dashed line). The upper panel shows the case when the mass 
accretion is limited by a maximum free fall time of 1000 seconds, while the lower panel shows 
the results for no limiting tff (cf. Fig. CD). 
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Fig. 7. — The value of critical specific angular momentum during an iteration, for one 
exemplary value of the initial normalization parameter, x = ZoAcrit — 0.05. The models of 
the distribution function f(r,0) are: C (solid squares) and D (open circles), and accretion 
scenarios are: a (solid line) and c (dashed line). 
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Fig. 8. — The total accreted mass as a function of the initial normalization of the specific 
angular momentum distribution. The models of the distribution function f{r)g{6) are: C 
(solid squares) and D (open circles), and accretion scenarios are: a (solid line) and c (dashed 
line). The upper panel shows the case when the mass accretion is limited by a maximum free 
fall time of 1000 seconds, the middle panel shows the case of the specific angular momentum 
cut off to / max = 10 17 cm 2 s _1 , while the bottom panel shows the case of no free fall time limit 
and no specific angular momentum cut off. Note, that the solid lines for models C and D 
overlap. 
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Fig. 9. — The mass accreted with / > Z cr it as a function of the initial normalization of the 
specific angular momentum distribution. The models of the distribution function f{r)g{9) 
are: C (solid squares) and D (open circles), and accretion scenarios are: a (solid line) and 
c (dashed line). The upper panel shows the case when the mass accretion is limited by a 
maximum free fall time of 1000 seconds, the middle panel shows the case of the specific 
angular momentum cut off to Z max = 10 17 cm 2 s _1 , while the bottom panel shows the case of 
no free fall time limit and no specific angular momentum cut off. 
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Fig. 10. — The instantaneous mass accretion rate during an iteration, for the 4 models and 
2 chosen values of the initial normalization of the specific angular momentum distribution. 
The lower panel shows the models of the distribution function f(r): model A (solid squares) 
and model B (open circles) and the upper panel shows the models of the distribution function 
f(r)g(9): C (solid squares) and D (open circles). The accretion scenarios are: a (solid line) 
and c (dashed line). 



